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Since their discovery, Ultraluminous X-ray sources (ULXs) have attracted attention due to their combination of extreme 
luminosities and extra-nuclear locations. However, they are a fairly rare phenomenon, and attempts to investigate the 
general properties of the population have been hindered by a relative lack of known sources. Here, we present a large 
catalogue of ULX candidates including 655 detections of 475 discrete sources, based on the 2XMM Serendipitous Survey. 
To demonstrate the potential of such a resource, we present some scientific analysis of this population, focusing on the 
spectral turnover seen, often at ~6 keV, in the highest quality ULX data. We also demonstrate how the recent reflection 
and Comptonisation interpretations of this feature may be distinguished observationally in the future, specifically using 
NGC4517 ULX1, a previously unanalysed source with high quality data, as an example case. 
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1 Introduction 

Ultraluminous X-ray sources (ULXs) display the unusual 
combination of intense X-ray luminosities (exceeding 10 39 
erg s _1 ) and extra-nuclear locations, and as such have at- 
tracted considerable attention as attempts are made to iden- 
tify the cause of their apparently extreme emission. It is well 
accepted that the majority of ULXs are likely to be accreting 
black holes, and most explanations put forward are based 
on one (or more) of the following ideas: ULXs may contain 
larger black holes than stellar-mass black holes hosted by 
typical X-ray binaries (XRBs), possibly intermediate-mass 
black holes (Colbert & Mushotzky 1999; dynamical friction 
arguments almost certainly rule out the majority of ULXs 
hosting supermassive black holes; Miller & Colbert 2004), 
they may be stellar-mass black holes radiating above their 
Eddington limit (Begelman 2002), or they may be sources 
that do not emit isotropic ally, hence luminosities calculated 
assuming isotropy overpredict their power output (King et 
al. 2001). In the latter case, the optical nebulae observed 
around many ULXs seem to rule out extreme anisotropy 
due to relativistic jets orientated towards us (e.g. a 'micro- 
blazar'; Pakull & Grise 2008), but modest collimation by 
e.g. 'slim' accretion discs remains an intriguing possibility. 

Observationally, it has proven difficult to distinguish be- 
tween these explanations, notably between the IMBH and 
stellar mass interpretations. For a recent review of this topic, 
see Roberts (2007). In large part, this is because reliable dy- 
namical mass estimates are not yet available for any ULXs. 
However, ULXs are a fairly rare phenomenon, and a rel- 
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ative lack of known sources on which to base population 
studies has also hindered progress. Here, we present a large 
catalogue of X-ray detections of ULX candidates, as well as 
some analysis of its population, as a step towards addressing 
this issue. 



2 The Catalogue 

The main resource used in this work is the 2XMM Serendip- 
itous Survey (2XMM; Watson et al. 2009). In the interest of 
brevity, here we list only the main steps taken to condense 
2XMM into a catalogue of ULX candidates. For a full de- 
scription of the data reduction see Walton et al. (in prep.). 

- 2XMM was cross-correlated with the Third Reference 
Catalogue of Bright Galaxies (RC3; de Vaucouleurs et 
al. 1991), in order to select the X-ray sources within the 
elliptical D25 galaxy regions. 

- X-ray luminosities were calculated using the full band 
(0.2-12.0 keV) EPIC flux and galactic distances taken 
from Tully (1988) for galaxies with a recession velocity 
cz < 1000 km s _1 , or calculated assuming the Hubble 
flow otherwise (adopting Hq = 75 km s _1 Mpc" 1 for 
consistency between the two regimes). 

- Sources flagged as extended, and those fainter than 10 39 
erg s _1 by more than la, were discarded. 

- A minimum detection significance of 3.5er was required, 
to reduce the contribution from spurious detections. 

- Sources within 7.5" of the galactic centre were excluded 
to minimize the contribution of low luminosity AGN. 

- Known contaminants (e.g. AGN, SNe) were removed. 
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Fig. 1 Number distributions for detections in the presented catalogue, and some of its subsets, of the na deviation of the 
observed EPIC-pn band 5 counts from the expected counts calculated by extrapolating the lower energy spectral shape, 
assuming a simple absorbed power-law model. Left panel: distributions for the whole catalogue (top row), and the spiral 
and elliptical galaxy subsets (middle and bottom rows respectively). Right panel: distributions for the high, medium and 
low quality datasets (top, middle and bottom rows respectively; determined by band 5 EPIC counts, E5). Both panels: 
distributions in the left and right columns are calculated assuming (Nh) = 1 and 3 x 10 21 atom cm~ 2 . Most distributions 
are significantly skewed and extended towards negative values of n, which suggests the presence of a spectral turnover 
similar to that seen in the highest quality ULX data in a significant fraction of our population. 



These filtering measures produced a catalogue of 655 X- 
ray detections of 475 ULX candidates; 59 of these sources 
are nominally fainter than 10 39 erg s _1 but agree with this 
limit within their la luminosity uncertainties. Of these 475 
candidates, 307 are located in 142 spiral galaxies, with the 
remaining 168 located in 98 elliptical galaxies, despite the 
observations contributing to 2XMM covering twice as many 
ellipticals than spirals. This is most likely due to the link be- 
tween ULXs and star formation (Swartz, Tennant & Soria 
2009), as spiral galaxies typically display higher star for- 
mation rates than elliptical galaxies. 

Based on the work of Moretti et al. (2003), we have 
used the 2XMM sensitivity maps to estimate the contribu- 
tion from unidentified contaminants resolved from the cos- 
mic X-ray background (e.g. background AGN) within the 
catalogue, and find it to be only ~18 per cent. The frac- 
tional contamination is much higher for elliptical galaxies 
than spiral galaxies, due to a combination of them being 
larger in general, and so cover a greater cumulative sky area, 
and hosting an intrinsically less numerous ULX population. 

3 High Energy Turnover 

One of the key observational characteristics of ULX spectra, 
seen in the highest quality observations of ULXs to date, is 



that they display curvature in their high energy continuum 
(typically at ^5-6 keV); see Stobbart, Roberts & Wilms 
(2006) and Gladstone, Roberts & Done (2009). Such cur- 
vature is not frequently seen in the spectra of typical XRBs, 
hence may provide information on any physical differences 
between these types of X-ray source. In this section we at- 
tempt to investigate how widespread this spectral downturn 
is in the presented ULX population. 

A full spectral analysis of each detection included in 
the catalogue is beyond the scope of this work, so we turn 
instead to count rate ratios of the pre-defined 2XMM en- 
ergy bands, specifically between bands 3, 4 and 5 which 
cover 1.0-2.0, 2.0^1.5 and 4.5-12.0 kcV respectively. We 
adopt a simple absorbed power-law model, and generate 
empirical relations between the count rate ratios of bands 
3/4 and bands 4/5 and the photon index T with PIMMS 1 . 
These are used to estimate the band 5 count rate that should 
be observed assuming no curvature (constant T) from the 
band 3/4 ratio, which are in turn compared to the observed 
rates. An na agreement between the two is determined such 
that negative values of n imply fewer counts were observed 
than would be expected, i.e. that the spectrum does dis- 
play a downturn. We repeat our analysis with estimates for 
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Table 1 Values for the curvature parameter C (equation 
[TJ for the EPIC PN number distributions of n for the vari- 
ous catalogue subsets shown in Fig.Q] Negative values of C 
imply a spectral downturn is observed in the highest 2XMM 
energy band within the sampled population. 



Curvature Parameter, C (x 10 2 ) 



(Nh)" 


1 x 10 21 


3 x 10 21 


Whole catalogue 


-17.4 ±2.1 


-8.3 ± 1.7 


Spiral galaxies 


-21.4 ±2.9 


-10.1 ±2.3 


Elliptical galaxies 


-8.2 ±2.5 


-4.1 ± 1.7 


High quality 


-33.5 ±5.5 


-14.3 ±4.3 


Medium quality 


-7.5 ± 2.2 


-2.5 ± 1.3 


Low quality 


-11.2 ±2.8 


-8.1 ± 2.3 



a Column densities are quoted in atom cm 

the average neutral column density of (Nh) = 1 and 3 x 
10 21 atom cm -2 ; previous work suggests (Nh } should be in 
this range (e.g. Stobbart et al. 2006; Gladstone et al. 2009; 
Swartz et al. 2004) 

The EPIC-pn number distributions of n for the whole 
catalogue and some of its various subsets are shown in Fig. 
Q] We do not show the distributions for the EPIC-MOS de- 
tectors for clarity, but note that they display similar trends. 
The top row of the left panel shows the distributions for the 
whole catalogue, and it is clear that they are not symmet- 
rically distributed around zero, as would be expected if no 
particular spectral curvature was common in the population. 
Instead they are skewed and extended towards negative val- 
ues of 7i, implying the spectral curvature observed in the 
best quality ULX datasets may be seen in a significant frac- 
tion of our ULX population. In an attempt to quantify this, 
we define the curvature parameter C in equation Q] (where 
N is the number of detections); C should be consistent with 
zero for symmetric distributions. Values are quoted in Ta- 
ble Q]for all the distributions shown in Fig.Q] The presence 
of background AGN will have diluted these values, as they 
do not typically show curved spectra. In addition to the dis- 
tributions for the whole catalogue, the left panel of Fig. Q] 
shows a comparison between sources located in spiral and 
elliptical galaxies, while the right hand panel shows the evo- 
lution of these distributions with data quality, based on the 
number of band 5 EPIC counts. 

c= N(n>3)-N(n<-3) 
Ntot 

It is clear from both the right hand panel in Fig. Q] and 
Table Qjthat the proportion of the population which displays 
observable curvature increases with increasing data quality. 
This is as would be expected were the curvature intrinsi- 
cally present in the majority of ULXs; as shown in Stobbart 
et al. (2006) it can be subtle to detect even in the best qual- 
ity data. Many observations contributing to 2XMM would 
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Fig. 2 The EPIC-pn (black) 2-10kcV spectra of 
NGC4517 ULX1, modelled with an absorbed power-law 
(solid line; see text). This model clearly predicts an ex- 
cess of counts over those observed at high energies (above 
~6 keV). A data/model ratio plot is also shown in the lower 
panel. 

not have been long enough to reliably detect continuum cur- 
vature at high energies. The best quality subset of the cata- 
logue, which should provide the most reliable estimate, sug- 
gests this curvature is a common property of ULX spectra. 
At first, when comparing the populations from spiral and 
elliptical galaxies it may seem that high energy curvature is 
more common in ULXs located in spiral galaxies. However, 
it is likely that (Nh) is lower for the elliptical population 
than the spiral population. Taking this into consideration, 
we argue that high energy curvature is likely to be common 
to ULXs from both galaxy types. If so, given that ULX pop- 
ulations in spiral and elliptical galaxies are expected to be 
dominated by high- and low-mass binary systems respec- 
tively, the likelihood is that the curvature must be associ- 
ated with processes in the inner regions of the X-ray source, 
rather than any environmental effect. 



4 NGC 4517 ULX1 

One source that came to light during the production of the 
catalogue, located in the edge-on spiral galaxy NGC 4517, 
is NGC 4517 ULX1. The archived EPIC spectrum contains 
^15500 counts, so the data meet the quality criterion of Sto- 
bbart et al. (2006) and can be considered amongst the best 
quality ULX data available. Following their example, we 
checked the 2.0-10.0 keV spectrum for curvature, compar- 
ing single and broken power-law models. Even allowing for 
absorption by a substantial column density of Nh ^8 x 10 21 
atom cm~ 2 , as suggested from inspection of the full 0.3- 
10.0 keV spectrum, the broken power-law model is strongly 
favoured, with a break energy of ^5.5 keV, implying the 
high energy continuum does indeed display curvature simi- 
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Fig. 3 A comparison of the Comptonisation and disc re- 
flection interpretations for ULX spectra; both have been ex- 
tended beyond the XMM-Newton bandpass up to 100 keV. 
The models agree well below ~10 keV, but there is a clear 
discrepancy of many orders of magnitude between them at 
higher energies, due to the presence of the Compton hump 
in reflection spectra. 

lar to that seen in the other high quality ULX spectra. This 
is demonstrated in Fig. [2] 

A number of explanations of this feature have been pro- 
posed, including Comptonisation of seed photons from the 
accretion disc in an optically thick corona, which shrouds 
the inner regions of the disc (Gladstone et al. 2009) and a 
combination of relativistically blurred iron emission and ab- 
sorption features, originating from the inner regions of the 
disc (Caballero-Garcfa & Fabian 2010). Other explanations 
invoke advection of radiation in the inner regions of a slim 
disc (Watarai, Mizuno & Mineshige 2001), although these 
have not been as successful at modelling the broad band X- 
ray spectra (see Gladstone et al. 2009). Both the Compton- 
isation and disc reflection interpretations are applied to the 
0.3-10.0 keV spectrum and found to provide statistically 
equivalent representations of the data. We do not provide 
full details of the spectral modelling here (see Walton et 
al. 2010, submitted), but merely comment that the results in 
each case are similar to those obtained from the application 
of these models to other ULXs: the electron temperature of 
the corona is low and its optical depth is high if Compton- 
isation is correct, while the iron abundance of the disc is 
highly super-solar and the corona very compact (within a 
few Rq) if reflection is correct. 

A detailed physical consideration of these results is pre- 
sented in Walton et al. (2010), but neither model may cur- 
rently be excluded, so other observational means are re- 
quired to distinguish between them. In Fig.[3]the two mod- 
els are compared after being extended beyond the XMM- 
Newton bandpass to 100 keV. Above 10 keV the Comp- 
tonisation model continues to curve downwards, while the 
reflection model turns back up. This is due to the presence 
of the 'Compton hump' in reflection spectra, a broad emis- 



sion feature often observed at ~30 keV, which arises due 
to the interplay between photoelectric absorption of low en- 
ergy photons and Compton down-scattering of high energy 
photons within the reflecting medium. The presence of this 
feature leads to a large difference in the expected emission 
predicted above ^10 keV by the two models. Currently the 
Comptonisation model is calculated assuming a purely ther- 
mal distribution of electrons, although it is likely the corona 
also contains a non-thermal population. Eventually Comp- 
tonisation from this population will arrest the downward 
curvature, however this may not become significant until en- 
ergies > 100 keV, and we still expect there would be large 
differences between the two model predictions. Hence, even 
fairly short observations of ULXs above ~10 keV could be 
very important in determining the underlying physical pro- 
cesses behind the observed emission from ULXs. 

With current instrumentation such observations are only 
possible for M 82 X-l due to the lack of high resolution 
X-ray imaging spectrometers that operate at such high en- 
ergies. The indication from M 82 X-l is that the spectrum 
continues to curve downwards, which favours the Compton- 
isation interpretation for this source (Miyawaki et al. 2009). 
The launch of the hard X-ray imagers aboard NuSTAR 2 and 
Astro-H 3 over the next few years should enable such obser- 
vations to be carried out for a large number of ULXs. 
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